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We explore the possibility to excite the so-called embedded solitons in specially designed slot waveguides 
based on silicon and silica or silicon nanocrystals. This requires the excitation of the structure with quasi-TM 
polarized pulses - for which Raman effect is absent - and at a specific infrared wavelength for which only the 
second- and fourth-order group velocity coefficients are non-vanishing. Pulses launched in these conditions will 
generate a spectrally localized continuum, associated to a frequency comb coming from the spectral interference 
of many embedded solitons. © 2012 Optical Society of America 
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In 1993, Hook and Karlsson [1] demonstrated theoret- 
ically the existence of a family of solitary waves, solution 
of the equation 

id z A - + ^!3 4 dfA + 7 \A\ 2 A = 0, (1) 

where A(z,t) is the electric field envelope, ^2,4 are re- 
spectively the second- and fourth-order group velocity 
dispersion (GVD) coefficients of the waveguide calcu- 
lated at a reference frequency uo$ (taken to be the central 
frequency of the input pulse), and 7 is the nonlinear co- 
efficient of the fundamental fiber mode. Eq. (1) is valid 
when ujq is located at a local maximum (or minimum) of 
the GVD, where /3s vanished identically. The properties 
of such peculiar solitary waves - sometimes called em- 
bedded solitons (ESs) due to their coexistence with linear 
waves [2] - have been since then studied in detail. It was 
found that ESs exist only for ^2,4 < 0, can have or not 
oscillating decaying tails, can form symmetric (in-phase) 
and antisymmetric (out-of-phase) bound states, and that 
only the single peak solitons can be stable [3-6]. 

In principle, it would be possible to design specific 
solid-core photonic crystal fibers (PCFs) [7] with a quar- 
tic profile of the GVD as required by Eq. (1), by using for 
instance the dispersion-flattened fibers reported in [8], 
or the nano-bore fibers reported in [9]. However, the 
perturbation induced by the Raman effect would com- 
pletely destroy the ESs. Indeed, from the very begin- 
ning of the propagation, their central frequency would 
continuously translate due to the intrapulse Raman self- 
frequency shift [10], thus moving all solitons away from 
the point at which ^3 = 0, which immediately breaks 
the validity of Eq. (1). This is mainly the reason why 
the above theoretical works have been considered only 
as a mathematical curiosity, only marginally discussed in 
textbooks [2,11], and largely forgotten during the past 
two decades. 

However, the recent advent of silicon photonics [12-15] 
has opened many interesting possibilities in nonlinear 
optics, mainly due to its potential applications in the 
spectral region extending from the near- to the mid- 
infrared. The high refractive index of silicon, combined 



with the silicon-on-insulator (SOI) technology, allows for 
a tight confinement of the optical modes and a conse- 
quent increase of the nonlinear coefficient of the waveg- 
uide, enabling efficient nonlinear optical interactions at 
low power levels and in relatively short length. 

In this Letter we propose (for the first time to our 
knowledge) a feasible way to observe and use the gen- 
eration of ESs in specially designed silicon-based waveg- 
uides. 

The two structures that we propose - both of which 
are built on a silica substrate - have been inspired by 
Refs. [16,17] and are depicted in Fig. l(a,b). In Fig. 1(a) 
we show the structural parameters of our silicon-silica 
(Si-Si02) waveguide, where the GVD for quasi-TE and 
-TM modes and mode profile for quasi-TM mode are 
also shown. In Fig. 1(b) we show the same information 
as in Fig. 1(a), but for our silicon-silicon nanocrystals 
(Si-SiNc) waveguide. The gray-shaded circles, centered 
at the wavelengths Ao at which Ps(Xo) = (which from 
now will be referred as quartic point), are the regions 
around which the GVD can be approximated by includ- 
ing only /3 2 and ^4 in the dispersion. Note that, in those 
regions, the dispersion is anomalous (#2 < 0) and the 
GVD curvature is negative (^4 < 0). As described above, 
these are exactly the necessary conditions for the possi- 
ble existence of localized ESs, and the parameters of the 
structure must be carefully designed in order to achieve 
the correct sign and curvature of the GVD at the desired 
wavelength. 

Apart from the shape of the GVD, one must consider 
the impact of the two-photon absorption (TPA) gener- 
ated by the free carriers around Ao- As it was demon- 
strated in numerous experimental works (see e.g. Refs. 
[18-21]), the impact of TPA is greatly reduced if A > 2.2 
/im, away from the two-photon band edge, where the 
three-photon absorption is also insignificant [19-21]. 

The final crucial property for the present work is that, 
for an SOI waveguide fabricated along the [110] direc- 
tion on the [110] x [001] surface, stimulated Raman scat- 
tering (SRS) cannot occur when an input pulse excites 
the quasi-TM mode of the waveguide [15,22]. Thus, for 
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Fig. 1. (Color online) The silicon-based structures (both on 
a Si02 substrate) studied in this Letter, (a) GVD of quasi- 
TE (dashed line) and -TM (solid line) modes of the Si-Si0 2 
slot waveguide. Inset below: TM-mode profile at Ao = 2592 
nm. Inset above: geometry and materials of the structure. 
Parameters are: h — 700 nm, w\ — 220 nm, W2 — 60 nm, 
w 3 = 600 nm. (b) GVD of quasi-TE (dashed line) and -TM 
(solid line) modes of the Si-SiNc slot waveguide. Inset below: 
TM-mode profile at Ao = 2511 nm. Inset above: geometry 
and materials of the structure. Parameters are: h — 700 nm, 
wi — 190 nm, u>2 — 80 nm, W3 — 570 nm. In both figures, the 
red dot indicates the wavelength (Ao) at which /3s = 0, and 
the gray area is the region where the quartic approximation 
of the GVD, and therefore the use of Eq. (2), is valid. 



pulses launched close to Ao = 27tc/cjo, the nonlinearity 
is dominated by the Kerr effect, and is not affected by 
the carrier dynamics, and the dimensionless propagation 
equation can be thus written as 

idrf + \d 2 T i> - |5 4 |^ + (1 + isd T ) M V = 0, (2) 



where we have performed the typical scalings A = v^A? 
t = tor, z = zo£, and where zq = to/l/^l, ^4 = 
^4/ (24 1/^2 1 tg), Po = (7^0) _1 ? 7 is the nonlinear coef- 
ficient of the structure calculated with a mode solver 
(mainly due to the field inside the central silicon part), 
and s = (cjo^o) -1 , with to being the input pulse dura- 
tion [11]. 

In Eq. (2), it was very important to take into account 
the shock operator that multiplies the Kerr nonlinear 
term. Such operator considerably influences the dynam- 
ics in those situations where the Raman effect is ab- 
sent, as in the present case. This constitutes the main 
difference between Eq. (1) and Eq. (2), the latter hav- 
ing now a direct physical interpretation in the systems 
studied here. It is therefore crucial to check that the 
shock term does not destroy the ESs around the quartic 
point. In order to verify this, we launched an input pulse 
^in = Nsech (r) in the simulation of Eq. (2), and propa- 
gate for several dispersion lengths, see Fig. 2, in absence 
[Fig. 2(a)] and in presence [Fig. 2(b)] of the shock op- 
erator, for the parameters indicated in the caption. It 
is clear by comparing these two figures that the shock 
term has the effect modifying the velocity of each indi- 
vidual ESs, without affecting their shape or stability dur- 
ing the propagation, and actually allowing the generation 



of more solitons than what occurs in the case of s — 0. In 
order to verify this fact, we compare the approximate an- 
alytical non-oscillating solutions of Eq. (2), for 8 = 0, ob- 
tained via variational methods, with the actual solitons 
produced in the simulations (with 8 ^ 0). Such solutions 
are given by [5]: ^(£,t) = f(r)e iq t, f = a(q)sech 2 [k(q)r}, 
with a(q) and k(q) being the two functions 
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The dashed line in Fig. 2(c) shows the relation between 
k(q) and a(q) given by Eqs. (3-4), while the dots indi- 
cate the parameters of the solitons extracted with a fit 
from the simulation given in Fig. 2(b) in the (a,fc) space. 
One can see that the dots are very well predicted by the 
dashed line, indicating the formation of true ESs in the 
simulations. 
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Fig. 2. (Color online) (a) Propagation of an input pulse ^in = 
iVsech (r), with N = 10 and S4 = —0.01, and no shock term 
(s = 0). Due to the non-integrability of Eq. (2), which induces 
the presence of a solitonic friction, the pulse splits into several 
fundamental embedded solitons, approximated by Eqs. (3- 
4). (b) Same as in (a), but in presence of shock term (s = 
0.03). The shock term induces an asymmetric splitting of 
the solitons, which all decrease their velocity with respect 
to the case of (a), (c) Position of the parameters (a,k) of 
each solitons formed in (b) after a propagation distance £ = 
10, on the predicted curve calculated by using Eqs. (3-4). 
It is evident that all the localized structures formed [which 
are numbered in (b,c,d)] are embedded solitons. (d) XFROG 
spectrogram of (b). 
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Finally, in Fig. 3 we show the propagation of a to = 40 
fs pulse in the structure of Fig. 1(a), for Ao = 2.6 /mi, 
with other parameters reported in the caption. The for- 
mation of a frequency comb extending from approxi- 
mately 2 to 3 jam is the consequence of the strong spec- 
tral interference between many ESs forming near the 
quartic point of the GVD, which are unable to sepa- 
rate in the frequency domain due to the absence of the 
Raman effect (and thus of any deceleration that can act 
on the solitons, followed by a continuous self-frequency 
shift, [10]) for the chosen quasi-TM input polarization. 
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Fig. 3. (Color online) Formation of an ESs-induced frequency 
comb. The bottom part shows the GVD (dashed-dotted line) 
of the quasi-TM mode vs wavelength for the structure shown 
in Fig. 1(a). The top part shows the input pulse spectrum 
(dotted line) and the output spectrum (solid red line) after 
a propagation of £ = 10. Input pulse duration is to = 40 
fs, input wavelength is Ao = 2.6 /xm, at which 02 — —0.05 
ps 2 /m and /?4 c± — 2.5 x 10 -5 ps 4 /m. The nonlinear coefficient 
for the structure of Fig. 1(a) is calculated by using the vector 
model of Ref. [23], and is 7 ~ 42 m _1 W _1 . The corresponding 
fundamental power is Po = 0.74 W. Input pulse shape is 
A/sech (t/to), with soliton order N = 5. 
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In conclusion, we have demonstrated analytically and 
numerically the possibility of the observation of localized 
embedded solitons in specially designed silicon-based slot 
waveguides. The absence of the Raman effect for TM- 
polarized pulses, together with specific structures that 
allow a quartic GVD curve to be realized far from the 
two-photon absorption region of silicon, permits the gen- 
eration of many non-radiative ESs and the formation of 
a spectrally localized frequency comb around the quartic 
point of the GVD. 
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